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D . Abstract 
^: 

It is shown that the mixing of lepton doublets of the Standard Model can 
yield sizeable contributions to the lepton asymmetry, that is generated through the 
^ . decays of right-handed neutrinos at finite temperature in the early Universe. When 

calculating the flavour-mixing correlations, we account for the effects of Yukawa as 
l/^ ■ well as of gauge interactions. We compare the freeze-out asymmetry from lepton- 

. doublet mixing to the standard contributions from the mixing and direct decays of 

\ right-handed neutrinos. The asymmetry from lepton mixing is considerably large 

when the mass ratio between the right-handed neutrinos is of order of a few, while 
it becomes Maxwell-suppressed for larger hierarchies. For an intermediate range 
between the case of degenerate right-handed neutrinos (resonant Leptogenesis) and 
the hierarchical case, lepton mixing can yield the main contribution to the lepton 
asymmetry. 



1 Introduction 

Sources for CP-violating effects are often categorised into contributions from mixing 
and from direct decays. This applies to Leptogenesis [ij as well, where usually the 
mixing (ij and the direct decays of right-handed singlet neutrinos iVj are accounted for. 
Of particular interest is the source from mixing, because in the situation where the mass- 
difference of the right-handed neutrinos is small, it gives rise to a resonantly enhanced 
CP- asymmetry (2- 10 



Besides the right-handed neutrinos, Higgs bosons and Standard Model lepton dou- 
blets directly take part in Leptogenesis. Direct decay asymmetries from Higgs bosons 
and lepton doublets are usually not considered, because in the vacuum, these particles 
cannot decay into heavy right-handed neutrinos. Consequently, at zero temperature. 



(A) (B) (C) 

Figure 1: Diagrammatic representation of contributions to the collision term that yield 
a source for the lepton asymmetry (A) through mixing of right-handed neutrinos, (B) 
through mixing of Higgs doublets and (C) through mixing of lepton doublets. Solid 
lines with arrows are propagators for the Standard Model lepton doublets £, dashed 
lines with arrows are propagators for Higgs doublets cf) and solid lines without arrows 
are propagators for the right-handed neutrinos A^. The double lines represent resummed 
propagators that in particular account for the mixing of flavour. 



there are no kinematic cuts through the propagators of a right-handed neutrino and a 
lepton- or Higgs-doublet that can lead to CP-violation. At finite temperature however, 



these cuts contribute to the asymmetry [llHl3j. since the plasma continuously absorbs 
and produces leptons and Higgs bosons. 

Concerning mixing, right-handed neutrinos, as these are sterile, are well suited for 
the efficient production of asymmetries since fiavour off-diagonal (where we define the 
fiavour of these particles by their mass-eigenstates, in contrast to Standard Model lep- 
tons) correlations are only damped by the Yukawa couplings, which are small by the 
requirement of a substantial deviation of the right-handed neutrinos from equilibrium. 

Gauge interactions, in particular in the Standard Model, are typically much larger 
than the Yukawa couplings of the right-handed neutrinos. However, as these are fiavour 
blind, they do not directly damp off-diagonal flavour correlations. Therefore, also the 
gauged particles within the Standard Model or its extensions can lead to substantial 



CP- violation from mixing. For multiple Higgs doublets, it is shown in Ref. [ij], that 
their mixing can be a viable source for Leptogenesis. In the present work, we show that 
the mixing of the lepton doublets of the Standard Model contributes a sizeable amount 
of CP-violation to Leptogenesis as well. 

In order to calculate the generation of the asymmetries in the flnite temperature 
background, we use the Closed-Time-Path (CTP) method js-lO, 13-31|. The diagram- 
matic representation of the source terms for the asymmetry in Figure [1] indicates how 
the mixing of the various species yields contributions to Leptogenesis . When we assume 
that the right-handed neutrinos are heavier than the lepton doublets and Higgs bosons, 
the imaginary parts of the self-energies that appear within the resummed propagators in 
Figure [U^B) and[T]^C) are purely thermal, as the corresponding 1^2 processes are kine- 
matically forbidden in the vacuum. This resembles the situation for Leptogenesis from 



the mixing of right-handed neutrinos with masses of the GeV scale or belo w |lOl . |32|-|35 



where the CP- violating cut is dominated by purely thermal contributions [lO|. For the 



direct CP-violation, the contribution of such non-standard cuts is calculated in Ref. 13 



In Section [21 we derive the off-diagonal correlations within the distributions of the 



2 



left-handed lepton doublets. Even though gauge interactions do not directly violate 
flavour and damp its correlations, they catalyse the decay of flavour correlations in the 
presence of flavour- dependent masses or interactions. For flavoured Leptogenesis, the 



relevance of such effects is discussed in Ref . [27[ . Here, we show that gauge interactions 



effectively limit the maximal resonant enhancement of a CP-asymmetry resulting from 
the mixing of gauged particles. For the mixing of Higgs doublets, this aspect is not 
discussed in Ref. IJ] , a shortcoming that is addressed in Section [51 

An important aspect within the present context is the role of lepton-number vio- 
lation and conservation. The spinor trace of diagram [T](c), that quantifles the rate of 
CP-violation from mixing of lepton doublets, contains no lepton-number violation that 
could be described by insertions of an odd number of Majorana mass-terms. Therefore, 
lepton number must in some sense be conserved. This, we show in Section [3l where it 
turns out that the source for the asymmetry in the lepton doublets cancels the helicity 
asymmetry within the right-handed neutrinos. The latter is however rapidly violated by 
the Majorana masses, provided the right-handed neutrinos are non-relativistic, as it is 
the case in the strong washout regime. 

In Section HI we then calculate the freeze-out value of the asymmetry that results 
from mixing lepton doublets in the strong washout regime. This, we compare to the 
asymmetry from the standard sources, i.e. from the mixing and direct decays of right- 
handed neutrinos. For illustrative purposes, we choose a particular slice in parameter 
space and vary the mass ratio of the right-handed neutrinos between one (degenerate 
case) and a few. It is a characteristic feature of the asymmetry from lepton mixing, that 
it is exponentially small for large mass ratios, as the CP- violating cut is purely thermal 
and becomes Maxwell suppressed in that situation. 

We close with summarising and concluding remarks, that are given in Section O 



2 Off-Diagonal Correlations of Lepton Doublets 



The model that we consider here is the same that underlies flavoured Leptogenesis [27 



36l438l| and is given by the Lagrangian 



C =^i^m{i^ - Mm^Nj + ^£ai#£a + {d^<p^){d^<P) (1) 

The Standard Model Higgs doublet is denoted by 0, the lepton doublets by £a, right- 
handed charged leptons by R^ and right-handed neutrinos by Ni. The associated spinors 
are denoted by ip with the subscript of the respective fleld. We use the deflnition 0^ = 
(pleba, but suppress explicit indices of the SU(2)l gauge group in the following. We 
transform the Ni to a basis where Mn is diagonal and write M^i = M^u. Furthermore, 
we take hab to be diagonal, where a,b = e, fi, r. We are mainly interested in the dynamics 
within the temperature range 10^ GeV < T < 10^^ GeV, such that only the r-lepton 
Yukawa couplings are in equilibrium. In that case, we can neglect h^e and /i^^ and then 
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perform a unitary transformation of the lepton flavours £e and to £o-± and such that 



y2(Tx = ^3ctx = 0) see Ref. 27| for the exphcit construction of such a transformation. 
(After this transformation, h is no longer diagonal within the block of the flavours e 
and r. As we may however neglect these Yukawa couplings at temperatures above 
10^ GeV, h can still be considered as effectively diagonal). Then, only lepton flavour 
asymmetries within and are produced, while there is no asymmetry in the remaining 
linear combination £a±, since iai. couples only to one of the right-handed neutrinos, and 
no CP- violating interference terms occur. While it is the salient feature of flavoured 
Leptogenesis that all CP-odd off-diagonal correlations between the flavours r and a 
are rapidly erased, because the interactions mediated by hrr are in equilibrium, we find 
in the present work that the deviation of the right-handed neutrinos from equilibrium 
does generically sustain CP-even [cf. Eq. fl22|) below] off-diagonal correlations of the 
leptons [Eq. f|T71) below], which in turn enter the source term for the CP- violating lepton 
asymmetry. We also note that the present discussion does not directly apply but can 
easily be generalised to temperatures below 10^ GeV, where the flavours e and /i are 
distinguishable. 

The present calculations based on the CTP formalism require numerous definitions 



271 ]. In particular iS*^'^ denotes fermionic and 



that can be found within Refs. 
iA*^'^ scalar Wightman functions. Subscripts on these indicate the particular fields and 
flavour correlations. Distribution functions for particles and anti-particles are denoted 
by / and /, number densities by n and n. Again, fields and flavour correlations are 
indicated by subscripts. A prefix 6 implies that the particular quantity is the deviation 
from the equilibrium propagator, distribution function or density. 

The dynamics of the out-of equilibrium distributions of left-handed lepton doublets 
can be described by the kinetic equation 27 



9hab 



-tr 



iC, + iC 



- ab 



(2) 



where rj is the conformal time in the expanding Universe. The functions g^^ " appear 



within the spinor components of the left-handed lepton propagator 



<,> 



into statistical and spectral terms in the narrow width limit [8|, [181, Il9|, |39|, |40 



and factorise 
When 



Ikl > VGT, where C 



m)9i - 



{l/2)gl, we may approximate 
---27rS{e)\km^)Sfeabik), 



(3a) 
(3b) 



For |k| <C \fGT, this approximation, in particular the relation between the sign of 
and the helicity, becomes inaccurate due to the presence of hole modes. However, this 
region of soft momenta covers only a small part of the phase-space, such that the approx- 
imations ([3]) are sufficiently accurate for the present purposes. The flavour-dependent 
self-mass terms are given by 



<b{k) 



16|k| 



(4) 
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where we have neglected the contributions from the couphngs Y. 
It is useful to decompose the collision term as 

c, = cf'''' + cj + c'; + c', (5) 

into contributions mediated by the Yukawa couplings Y and h and the gauge couplings g. 
Within Cj, we collect the CP -conserving processes mediated by Y, which are to leading 
approximation quadratic in Y. The term Cf'^^ is mediated by Y as well and accounts 
for CP-violating processes. To leading approximation, it is fourth order in Y. 
The part of the lepton self-energy mediated by Y is 



(^) = '^a^y,b I -§^,P^^S%{p)iAlf [p - k) , (6) 

where /, g are CTP indices. The term Cj is of importance for two reasons: First, as in the 
standard scenario of Leptogenesis, its diagonal components describe the washout of the 
CP-odd lepton asymmetries accumulated in the diagonal correlations of lepton flavours. 
(For flavoured Leptogenesis, the off-diagonal correlations are rapidly erased due to the 
Yukawa coupling hrr-) Second, its off-diagonal components correspond to the source of 
CP-even, flavour off-diagonal correlations in the lepton doublets. As the CP-odd off- 
diagonal correlations are suppressed by the Yukawa interactions hrr, and the CP-even 
off-diagonal correlations are suppressed by an order of Y'^ times the deviation of the 
right-handed neutrinos from equilibrium, when compared to the diagonal distribution 
functions, we can substitute the flavour diagonal components of 15"/'^ for the purpose of 
calculating the contributions to Cj that are of leading importance. We thus obtain 

CLik)+Cfr = it]aAk)iS<,{k) - iti^Ak)iSl,{k) (7) 
^ - [Y^YUPn I ^ [iS>,,{p)iA<{p - k)iSUk) - iS<,,{p)iA>{p - k)iS>,{k)] 

j_nCPV 

where we take iA^'^ to be of equilibrium form. Note that while the off-diagonal corre- 
lations within i and A^, which are always out-of-equilibrium, can be neglected for calcu- 
lating Cj to second order in Y, within the CP violating source C^^^, these off-diagonal 
correlations are essential for the asymmetry that we calculate in Eq. fl27|l below. 

For definiteness, we assume strong washout, i.e. Mj^i ^ T at the time of Leptoge- 
nesis. (Due to the exponential Maxwell-suppression, masses of a factor of a few above 
the temperature are enough for sufficiently accurate approximations.) This allows us to 
neglect quantum statistical blocking and enhancement factors on many occasions and 
to use Maxwell distributions instead of Fermi-Dirac or Bose-Einstein distributions. We 
parametrise the deviation of the neutrino distribution from equilibrium by a pseudo- 
chemical potential fi^i as 

Sfm{p) = ^e"^V^^.. (8) 
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While in the present setup, there are no interactions that directly force kinetic equi- 
librium of the sterile neutrinos, it turns out that above parametrisation indeed gives a 
good approximation to the distribution function of sterile neutrinos in strong washout 
scenarios 



26, 41 



We now aim to compute the off-diagonal distributions f^ab and fiab, a ^h. In kinetic 
equilibrium, these can be inferred from the number densities of leptons and anti-leptons. 



oo 









<5»7«. = -2 / ^/^SSW. (9b) 



-oo 



(Note that ± on 5?^-^^ refers to the particle and anti-particle, whereas on g'^'^, it refers 
to helicity.) The prefix 5 indicates that these quantities are deviations of the number 
densities from their equilibrium values. Consequently, the charge density q^ab = ^''^^ab ~ 
(5n^{, can deviate from zero. Pair creation and annihilation processes tend to strongly 
suppress the combination Sn^^^ + ^^iab- However, it turns out that for small Yukawa 
couplings, the finite value of this combination limits the maximal resonant enhancement 
of the asymmetry, cf. Eq. f ll7p below. For this reason, equations for 5n^^ rather than 



just qgab need to be derived, as a generalisation of the methods introduced in Ref. [14 . 

In order to obtain these densities when the time-derivative in the kinetic equation ([2]) 
is small compared to the effective mass difference Cfa ~ Cbb^ useful to multiply these 
equations by |k| and then to integrate over d^k. Doing so, we encounter on the right-hand 
side 

oo 

- / ^ / - / f / - - E i^.^««r . (10) 

-oo ' 

In the limit of strong washout, where the quantum statistical distributions can be ap- 
proximated by Maxwell distributions, we obtain 

^,^_rlMim-^^ (11) 

2 2 712 I 

Clearly, this term is exponentially suppressed for large M^i. It can be interpreted as 
the CP- violating cut that appears in the resummed lepton propagator of Figure [T]^C), 
and since the reaction £ -H- 0* -f iVj is kinematically forbidden in the vacuum, this cut is 
purely thermal. 

The interactions mediated by the Yukawa couplings h lead to the decay of off-diagonal 
correlations of the lepton doublets. The relevant contribution to the collision term is 
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given by [27 

CUk) ^hlK, I ^ [iS>Jp)iA>{k -p)- iS<M^A<{k - p)] i6S,,,{k) . (12) 
In the integrated kinetic equations, then the following terms occur: 

oo oo 

dk^ f d^k ,^ , . f dk^ f d^k f dS 





X 





tr [iS>M^A>{k -p)- iS<MiA<{k - p)] i6S,,,{k) 

oo 

, t , . + fdk' f d^k f d^p 



xtr [iS>Jp)iA>{k-p)-iS<M^A<{k-p)] 2St,{k) ^ 



=:-Bt[h^h5nj]^ 

and likewise 

° dfco f d^k 



— oo 

With these coefficients, the integrated kinetic equations are 

(15b) 

i 

In principle, there also occur terms 

<x{[Y^YU + [Y^YU)6ni, (16) 

on the right hand sides of these equations. However, we assume here that the dominant 
flavour-sensitive interaction is mediated by the lepton- Yukawa couplings. In particular, 
hrr dominates over the relevant elements of Y, which is valid provided Leptogenesis 
occurs at temperatures below 10^^ GeV. 
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The coefficient is the averaged rate for flavour- sensitive processes that are medi- 
ated by the couphng h. As the Higgs boson as well as the Standard Model leptons are 
massless in the symmetric Electroweak phase, these reactions require the radiation of 
an additional gauge boson. A systematic calculation is currently on the way, but from 



the similar problem of the production of light right-handed neutrinos [42|-|46[, we may 
estimate that sf = 1.0 x IQ-^T^, see the Appendix |X1 

Pair creation and annihilation processes may change Snf^^^. Even though these pro- 
cesses are flavour blind, we will see below that in interplay with flavour-sensitive pro- 
cesses, they can be important for the suppression of flavour off-diagonal correlations. 
Accounting for the large number of Standard Model degrees of freedom, these processes 
are dominated by the s-channel exchange of gauge bosons, and we estimate these within 
Appendix [A] as Bf = 1.7 x lO'^T^. 

When neglecting the time-derivative, Eqs. ( 1T5]) yield the induced off-diagonal (a 7^ b) 
correlations 

Qiab = qiabiv 00) (17a) 



■ [{hla - h'jTyiQf + (hi, + hl)Bf[2B^, + (hi, + hl)Bf] 
{Qeab/T')J2Yl,Y,>^Br, 

i 

2ihl + hl)BtZ^Y^^Y,,Br 



(17b) 

The general solution to Eqs. f|T5l) involves the damping rates 

rL = + {hi + hl)Bj ± ^JBf - [{hi - hDT^/lQf^ , (18) 

and is given by 

g^ab(r/) = qiab{ri ^ 00) + Qie-^i-6'' + Qae"^^'"''' . (19) 



This implies that it takes a finite time for these CP-even off-diagonal correlations to 
build up. In the present case, the rate for this is given by F"^^ ~ h^^T. This in 
particular implies that in the regime of flavoured Leptogenesis below 10^^ GeV, when the 
interactions mediated by h^-^- are in equilibrium, flavour off-diagonal lepton correlations 
are present, and their magnitude is given by Eqs. ( IT7I) . 

We express the resonant enhancement of the asymmetries through the factor 

Q ^ {hl-hl){T'/8) 

Waa - hDllQf + {hi + hl)B^[2B3 + {hi + hl)B^] ■ 
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The fact that 5nlab^^''hab suppressed but not vanishing imphes that kinetic equihbrium 
mediated by pair creation and annihilation processes is only approximately established. 
Processes of kinetic equilibration that do not mediate between 5rtl and 5nJ , i.e. which 
do not correspond to pair creation and annihilation, are however not conflicting with the 
solution ( |T71) . We therefore describe the distribution functions in terms of the number 
densities as 



Jk|/T 



T3 (elk|/T + l 



,2 ' 



I25n, 



2^3 



,|k|/T 



\2 ■ 



(21a) 
(21b) 



Notice also that from Eqs. 
n 



it follows that Ue 



"^labi ^lab 



n^* and therefore 



''tab 



77,^^^. This implies that the flavour off-diagonal lepton correlations are CP-even, 



c 



£ba\ 



p)] C^ + 0(Y^ 



(22) 



such that to this end, the lepton asymmetry is not yet broken. 



3 Asymmetries of Lepton Doublets and Singlet Neu- 
trinos 

In this Section, we calculate the CP-odd charge asymmetries within the diagonal compo- 
nents of the distribution functions of lepton doublets as well as the helicity asymmetries 
within the singlet neutrinos, that are induced by the off-diagonal, CP-even correlations 
of lepton doublets. In the strong washout scenario, helicity asymmetries within the right- 
handed neutrinos are quantitatively irrelevant, because of the Majorana masses, which 
effectively erase such asymmetries when the right-handed neutrinos decay. However, 



helicity asymmetries can be of crucial importance in weak washout scenarios [10|, I32|-|35 

In the present context, we are interested in the helicity asymmetries, because the 
lepton-doublet propagator and its self-energy correction do not encompass lepton-number 
violating combinations of operators, such that the total source for the lepton asymmetry 
plus the helicity asymmetry in the right-handed neutrinos should vanish. This constraint 
provides a useful consistency check that we perform in the present Section. In spite of 
the vanishing of the total source, the asymmetry of the singlet neutrinos is completely 
washed out in the decays, that proceed via the Majorana mass term within the strong 
washout regime. In contrast, the active leptons experience finite washout rates from the 
inverse decays, such that a total lepton asymmetry remains at freeze out. 
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The self-energy for the right-handed neutrino is js], [lol, 26 



(23) 



dp 



c 



where = 2 is the dimension of the SU(2)l representation of the lepton doublet. The 
self-energy enters into the collision term for the right-handed neutrino as 



trC^(A;) =tr it^^{k)iS<{k) - it'if{k)iS>{k) 



(24) 



A CP- and helicity-violating contribution C^^^ is induced by the flavour off-diagonal 
correlations of the lepton doublets. We denote this as the source term 



S 



Nij 



d^k 1 

(2^4*'' 



nCPy 



E 

a,b 



d^k 



{k) + cr\k) 

d^p 



(25) 



(27r)4 J (27r) 



i6Seat{p) {iA>{k - p) [iS<,,{k) + iS<^,{k)\ 
iA<{k-p) [iS>,,ik) + iS>^^ik)]) 



We have made use here of the (9(y^) lepton number conservation, Eq. f l22|) . and as- 
sumed that there is no charge asymmetry within the Higgs field, such that iA^'^(p) = 
iA^'^{—p). Moreover, the off-diagonal correlations within the right-handed neutrinos, 
which are important for the usually considered CP-violating sources for Leptogenesis, 
are neglected here. (An example with a quantitative comparison of the source from mix- 
ing and decays of right-handed neutrinos and the source from mixing lepton doublets is 
presented in Section HI) 

The source term can be expressed to enter the integrated kinetic equation for the 
right-handed neutrinos in the form 



-d„ 
2 ^ 



itT[^'S<'>{k)] = -S^,,. 



(26) 



The integral on the left hand side corresponds to the rate of producing a helicity asym- 
metry, which is in the limit of vanishing Majorana masses the same as an asymmetry 
between right-handed neutrinos and anti-neutrinos. Note that this integrated kinetic 
equation does not account for the total production of singlet neutrinos. In order to 
obtain these rates, the collision term is to be integrated from to 00 rather than from 
—00 to 00 over dk^ 26|. The factor | on the right hand side is a suitable normalisation. 



because in writing the Majorana field as a four component spinor, the number of degrees 
of freedom is doubled, though they are related by the Majorana constraint. 
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As for the lepton collision term Cj , the approximation ([7]) is used in order to calcu- 
late the CP-even off-diagonal correlations in terms of the deviation of the right-handed 
neutrinos from equilibrium. For that purpose, we have taken Sg to be of equilibrium 
form under the integral. Now, we want to obtain the CP-odd asymmetry in the lepton 
doublets that is induced by these off-diagonal correlations. It follows from integrating 
(Jtf as 



£ab 



d^k 1 



(27r)4 2 



tr 



C 



CPY 



E 

i,c 



{k)+Cf'"'\k) 
d'^p 



(27) 



(27r)4 J (27r) 



rtr 



(PRi5>,,(p)iA<(p -k)- PniS<M^A>{p - k)) 



X {Yl^,-i5S,,,{k) +i5Siac{k)YyY,, 



lb 



Besides this source, at fourth order in Y, there are also the standard contributions from 
direct CP- violation in the decays and inverse decays of right-handed neutrinos as well 
as from their mixing These add linearly to the asymmetry from lepton- doublet 
mixing, that is the main topic of the present work. For simplicity, we do not account 
for the standard sources of asymmetry in most of the following discussion and refer 
to Refs. [s], [lo|, 13, 26, 27] for the calculation of these contributions within the CTP 
formalism. In Section |4] however, we present a comparison of the asymmetry resulting 
from lepton mixing to the asymmetry that results from the standard sources. 
When expressing the kinetic equation for the lepton doublets as 



and observing that 



d^k 



^^^tT[Ce{k)+Ce{k)i 



tr[S, 



N 



-g.ujtr[Si 



(28) 



(29) 



we see that the sum of the helicity asymmetry within the right-handed neutrinos and 
the charge asymmetry in the lepton-doublets is conserved, unless taking proper account 
of the helicity-flipping Majorana masses M^. This is expected, since the self-energy 
that leads to the mixing of the lepton doublets, does not encompass lepton-number 
violating operators in the form of odd powers of the Majorana masses within the 
approximation represented by the diagram in Figure Wl^c) {i.e. when accounting for 
lepton-doublet mixing only and not for the mixing and the direct decays of the right- 
handed neutrinos). 

Now, since in the strong washout scenario, Mtvm ^ T, the right-handed neutrinos 
are non-relativistic, and they decay at tree-level with equal likelihood into leptons and 
anti-leptons. This is different for relativistic right-handed neutrinos, that are converted 
to leptons via scattering processes with the particles in the plasma. These reactions are 



11 



approximately helicity conserving. In the following, we calculate the asymmetry in the 
strong washout scenario, and for this reason, we only have to account for the production 
and the washout of the asymmetry within the leptons £, because the helictiy asymmetry 
within the right-handed neutrinos is rapidly violated by the Majorana masses. In the 
strong washout regime, all energies of particles that participate in decay and inverse 
decay processes are much larger than the temperature. One may therefore approximate 
Eqs. ( 12T|) by their Maxwell forms 



126nj 
I25n} 



|k|/T 



|k|/T 



(30a) 
(30b) 



We then substitute above approximations into Eqs. f[T7l) . and these into the source term 
for the lepton asymmetry (1271) . It enters the kinetic equations ((21), when integrating 
these over dk^. 



327r3 



^tv. 69fca Qqiac ^ 



T3 



Mf,^TK,{Mm/T) + 



(2vr)^ 



CLik) (31) 



167r3 T4 
d^k 



Here, we have evaluated the source term (1271) explicitly in the limit of strong washout, 
while the remaining integral corresponds to the washout term. 



4 Freeze-Out Asymmetries in the Strong Washout 
Regime 

We now calculate the freeze-out value of the asymmetry produced in the out-of-equilibrium 
decays of the individual iVj in the strong washout regime. The total asymmetry is 
then obtained as the sum of the particular asymmetries. As usual, it is convenient to 
parametrise the evolution by 

z, = Mm/T. (32) 
The scale factor a can be expressed through the comoving temperature as 

mpi 



a = aR?7 = Or 



Or 



Ni 



45 



(33) 
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where = 106.75 is the number of relativistic degrees of freedom of the Standard Model 
at high temperatures. The physical temperature is given by T = a^/a = Mj^ijzi. In 
order to take account of the expansion of the Universe, in all kinetic equations written 
down to this end, the mass terms should be multiplied by the scale factor a and all 
temperatures be replaced by the comoving temperature cr, c/. Ref 26|. The kinetic 
equations then take the form 



dzi 
dzi 



dYjsfk ^ eq^ 



(34a) 



(34b) 



where Yia = qeaa/s and Y^i = riNi/s, and s is the entropy density. The various distribu- 
tions, collision and source terms are (c/. Ref. l3 ]) 



9 ^"JVj 



3aR^/e M|^. [YY% 



223/27]-7/2 [YY^jj^^'^''^ ^aj^^cY^i^iaJ , (35a) 
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yeq _ ^-1/2 -3/2 
^ Nk — ^ ^ 



Nk 



M, 



Ni 



3/2 



^3/2 „3 „-2,Mjvfc/Mjv, 



/s 



^Nk - ^ YkaY^i^aRi 



M 



Nk 



(35b) 
(35c) 



Ni 
la 



5/2 



29/2^5/2 \ 



M 



Nk 



5/2 



E^TVi ^5/2^-^^Mjvfc/MA 



(35d) 



To the leading order _in^ deviations from equilibrium, Y^rj — Yj^- , the formal solution 



to Eqs. ([M]) is given by |47|, |48 



jJcNi dz' Ni 
'-'iaa 



rl-,l -' Ni •'. ea y ' 



c 



Ni 



(36) 



- f dz" Y.Biakz"^^'' cxp(-z"MNk/Mm) 

dz'a'^Hz')e ^' * 



where we define 



Ni dz' Ni 



I Ni 



c 



Ni 



M, 



Nj 



Ni 



(37) 
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For simplicity, we now consider Leptogenesis from two right-handed neutrinos A^jj. 
This is relevant in the situation when there are only two right-handed neutrinos in the 
theory, when additional right-handed neutrinos are much heavier than A^i and or 
when these decouple due to the smallness of their Yukawa couplings. 

The exponent in Eq. f l36l) is then extremal for z' = Zf with 



Mm + M, 



Nj 



Mm 



(38) 



which we solve numerically for Z{, what determines the freeze-out temperature. Laplace's 
steepest descent method then yields the flavoured freeze-out asymmetries (c/. Refs. j47|. 



m) 

YiTiz = 00) = aiaizf) 



2-K 



f Efe dz"Bfcfc2"'''' cxp{-^"Mjvfe/Mjvi) 



Finally, the total lepton asymmetry is the sum 



(39) 



(40) 



In order to assess the relative importance of the asymmetry from mixing lepton 
doublets, we compare it to the lepton asymmetry that results from the wave function 
and vertex corrections to the decays of right-handed neutrinos TV. We can simply express 



dz 



(41) 



where again, Y/^^ is the asymmetry that results form the decay of Ni. The decay asym- 
metry is 0, 49 



where Xj = M^j/M^i and 



Im 



yW-iYW- 



3(x, - 1: 



e(x) 



-X 



'I + x) log 



1 + X 2 — X 



X 



X 



The freeze-out temperature is now determined by 



(42) 



(43) 



(44) 
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where Bfjf, is given by Eq. (I35dp . The factor that appears in the formal solution (l36l) 
for y*^* is now 



af (z,) = -ef ^2-^7r-^z?e-^^ . (45) 

We now pick some particular, but illustrative points in parameter space in order to 
compare the asymmetry from the mixing of lepton doublets to the asymmetry from the 
wave function and vertex corrections to the decay parameter of right-handed neutrinos. 
Parameters consistent with neutrino oscillations can be constructed with the help of the 
relation 5^ 



3 



3 



(46) 



^ V 

where m^'^^^ = diag(mi, m2, ms) is the diagonal mass matrix of active leptons, = 
diag(M7vii, M^22, M^ss) and v = 246 GeV. The PMNS matrix is 

with U±s = diag(e^*''/^, 1, e"*^**^/^). The non-zero entries of the matrices \/(*-') are 

V^^) = Vff^ = cos9,,, = sin%, V^f = -sin%, V^^ =,1. (48) 

The 9ij are the mixing angles of the active neutrinos, and ai, a2 and 5 are CP- violating 
phases. The matrix TZ must satisfy TV'TZ = 1. We parametrise it by the complex angles 
Uij and the matrices TZ^^^^ with non-zero entries 

= cosc.,,, 7^^f = sinc.,,, 7^(f = -sinc.,,, 7^i|\= ^ 1 , (49) 

Tl = 7^(23)^(13)^(12) _ 

We take Mm3 ^ Af7vi.2 or alternatively, make the assumption that there is no right- 
handed neutrino N3 or simply that the Yukawa couplings of N3 are negligibly small. The 
parameter space can then be restricted by the choice = 7r/2 and ^23 = 0. Additional 
parameters that we choose are listed in Table [TJ Finally, is a unitary matrix that 
brings y to a form, such that only asymmetries in two flavours r and a are produced, 



where a is a linear combination of e and /i, see Ref. [27[ for the explicit construction. 
(Explicitly, this means that Y2aj^ = Y^u^ = Y^a = Y^^ = 0, while the remaining entries 
are generically non-zero, and a± is the linear flavour combination orthogonal to r and 
a.) 

The one remaining parameter is M]y2, that we vary for values larger than M^vi- In 
Figure [21 we compare the resulting asymmetry from lepton-doublet mixing to the asym- 
metry from mixing and direct decays of right-handed neutrinos, where we take account 
of the decays of both, Ni and N2. For Miy2 M^i (degenerate regime), the asymme- 
try from the mixing of right-handed neutrinos is strongly enhanced, which is known as 
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Figure 2: Lepton asymmetry from mixing leptons (solid line) and from wave-function 
and vertex corrections for the right-handed neutrinos (dashed line). The model param- 
eters are listed in Table [TJ 



mi 





m2 


8.7 meV 


sin 9i2 


0.55 


sin 92'i 


0.63 


sin 6*13 


0.16 


5 






«! 





02 


Tx/2 


^12 


-7r/2 + i/2 


Wl3 


7r/2 


1^23 





Mni 


10^° GeV 



Table 1: Parameters for the numerical example in Figure El 



resonant Leptogenesis. The asymmetry from lepton-doublet mixing also is strongest for 
M]s[2 ~ M^i. For larger values of Mjv2, we clearly see the exponential suppression of 
the asymmetry, because the kinematic cuts that give rise to the CP-violation are purely 
thermal and are Maxwell-suppressed for large values of M]^2 (hierarchical regime), c/. 
Eq. (fTT]) . For intermediate values of M7V2, between the degenerate and the hierarchi- 
cal regimes, the asymmetry from mixing lepton doublets dominates over the usually 
considered asymmetry from mixing and direct decays of right-handed neutrinos. 

This example should illustrate that for M^ri ~ M]^2 (more precisely, when the right- 
handed neutrino masses are within a factor of a few), there generically is a sizeable 
contribution to the lepton asymmetry, that originates from lepton mixing. While we 
have considered here a particular slice of the parameter space, it would be interesting to 
perform a more systematic parametric study of Leptogenesis from mixing lepton doublets 
in the future. 
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Maximal Resonant Enhancement in Leptogenesis 
from Multiple Higgs Doublets 



In Ref. [IJ], Leptogenesis from mixing Higgs doublets is introduced. The model is based 
on the Lagrangian 

k 

- J2 Kki^l^i - J^^Ymk^NAPL^im + h.c.) . (50) 

kl mk 

Notice that the Yukawa couplings Y here are different from those in the Lagrangian ([1]). 
Using the same notations and conventions as in Ref. IJ], we obtain the off-diagonal 



correlations in the number densities of liiggs and anti-Higgs particles as 

oo,— oo 

Sntu = 2 J ^k'iDMk), (51) 



where is the Wightman function for the Higgs field. 

The kinetic equations for the off-diagonal correlations of the Higgs particles are then 
(we choose a basis where M^, that may include thermal corrections, is diagonal) 

±i(M|ii - Ml,,)6n%, = -[Y^Y],,Bl - J^iVn + vDH^i ' % (K12 + ^^^12) • 

(52) 

We have assumed here that the flavour-sensitive interactions are dominated by couplings 
Dik of the Higgs-fields to Standard Model fermions. Note however that the couplings 
Y and the quartic interactions among the Higgs fields are flavour sensitive as well and 
should be taken into account when the are not dominating. An estimate for the 
weighted annihilation rate 5^ 1.0 x 10~^T^ is given by Eq. ( l62l) in Appendix |Al As for 

-B^, we do not provide an estimate here, since it is model dependent. If one of the Higgs 
fields is the Standard Model Higgs field, this rate will be dominated by the top-quark 
Yukawa couplings yu. Because the scattering processes then also involve gluons, we may 
expect this rate to be enhanced by a factor of a few when compared to B'^. 
Eq. fl52|) has the solution 

g<Ai2 =^nl^2 - '^^<*i2 = 7^02i[r"^r]l25J , (53a) 

(53b) 



Vis ""-012 
^ _ ^011-^122 
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The enhancement factor IZ^ is analogous to the quantity fl20l) Q^ab for mixing leptons, 
but we have choose a different dimensionahty for convenience and the sake of comparison 
with Ref. Q. 



FoHowing Ref. we then obtain for the freeze-out asymmetry 

oo 



YpAz = oo) 



lm[Y,,Y;,Y,,Y*] 135 M^, 



= [Y^Y] 



tr[YYt] Uix^g^ 
9 / 5 mpi 



327r4 V 2^^ Mtv ' 



-\W., ((-2/5) X (2/5,)^/^) . 



(54a) 

(54b) 
(54c) 



(Note a missing factor 1/a/S^ in the expression for Ym in Ref. 14|.) 



6 Summary and Discussion 

We have calculated the mixing of Standard Model lepton doublets in the early Universe 
and studied its impact on Baryogenesis via Leptogenesis. The mixing is driven by non- 
equilibrium right-handed neutrinos, as described in Section |2l With no deviation from 
equilibrium, there is no mixing of the lepton doublets. The amount of mixing is quantified 
by the the correlations of the off-diagonal number densities (fT7|) . 

Substituting these off-diagonal flavour correlations back into the collision term for 
the leptons, we obtain the source for the lepton asymmetry, Eqs. ( 127|3T|) in Section |3l 
This asymmetry is opposite to the source for the helicity asymmetry in the right-handed 
neutrinos (which is however rapidly violated through the Majorana masses). 

By calculating the freeze-out value of the lepton asymmetry, we compare the source 
from lepton mixing to the standard source from mixing and direct decays of right-handed 
neutrinos in Section HI The standard result is enhanced for mass-degenerate right-handed 
neutrinos (resonant Leptogenesis), while the asymmetry from lepton mixing is Maxwell- 
suppressed for larger mass ratios. In between these regimes, we find that lepton-doublet 
mixing may be the dominant source of CP- asymmetry for Leptogenesis. 

Results for the impact of gauge interactions on the mixing of Higgs-doublets and im- 
plications for the freeze-out value of the induced asymmetries are presented in Section [51 

In the future, it would be interesting to perform a more systematic scan of parameter 
space, as the asymmetries presented in Figure [2] are from a particular slice only. For 
that purpose, a more systematic and accurate calculation of the rates of right-handed 
neutrino production and pair creation and annihilation of doublet leptons, that we es- 
timate here in Appendix |Al would be beneficial. The relevant techniques are currently 



under development, see Refs. 42|-|46|. 



The Closed-Time-Path method so far has already proven useful in order to describe 
the emergence of the baryon asymmetry of the Universe in a more systematic manner 
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and in order to calculate certain finite-ternperature corrections to the app roximations 



that are most commonly applied js- 10, 13, 20-29] . Together with Ref. 14], the present 
work shows that these methods can also be useful in order to identify sources for the 
CP-asymmetry, that have not been considered before. The further development of non- 
equilibrium field theory techniques and their application to the comparably simple prob- 
lem of Leptogenesis will therefore prove helpful for progress on both, the qualitative and 
quantitative understanding of the emergence of the baryon-asymmetric Universe. 
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A Pair Creation and Annihilation Rates and Flavour 
Sensitive Rates 



Pair creation and annihilation can be mediated either through the s-channel or the t- 
channel. Individual s-channel processes are typically suppressed compared to t-channel 
processes by a factor of a few (essentially the logarithm of the squared gauge coupling), 
because the latter may be infrared enhanced due to the exchange of soft particles. 
Nonetheless, due to the large number of matter degrees of freedom in the Standard 
Model, it should be a sufficient approximation to consider s-channel processes for pair 
creation and annihilation of matter particles and Higgs bosons only. Technically, it 
should also be possible to account for the t-channel contributions, which would however 
be more complicated due to the infrared enhancement and its screening. 

A standard procedure for calculating cosmological interaction rates of non-relativistic 
particles relies on the reduced cross section (t(s), where s is the Mandelstam s- variable. 
The main simplification arises due to the use of Maxwell statistics. For massless particles, 
this approximation incurs an order one inaccuracy. Moreover, the t-channel contribu- 
tions may leave the reduced cross section infrared divergent, a problem that should be 
accounted for by including the screening effects in the plasma. In conjunction with the 
comments above, we should expect that the following procedure yields an 0{1) estimate 
of the weighted pair creation and annihilation rate Bf. A more accurate calculation of 
this quantity may be subject of a future study. 

The total annihilation rate of massless particles A into massless matter or Higgs 
particles B in Maxwell approximation is given by 

^^^-^^-y 2|k^|2|k^|2|k^|2|k5|' ' \Maa^bb\ (55) 

00 
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The squared matrix element is understood to be a sum over all external polarisation and 
weak isospin states. The reduced cross section is defined as 



a(s) = — / dt\MAA^BB\- (56) 
Svrs J 



Let fA,B denote chiral fermions, bosons. We obtain that 

IfAfA^fBfB = , (57a) 

^fAfA-^bBbB = *^Y927r5 ' ^^^^"^ 

IbAbA^bBbB = ^384^ • (^'^^) 

When fx,Y and bx,Y are SU(2)l doublets, G = + gfilAlB)^ where qA,B denotes the 
U(l)y weak hypercharge of the particles. For singlets, G = gfiqAlB)- For the doublets, 
we have to divide the rates by two in order to average over the weak isospin such that we 
obtain the production rate for a single isospin state. Summing over all Standard Model 
degrees of freedom for B, we obtain 

_ / 97 4 163 A 4 

7«--.anything " [-^92 + -^^9l j ^ ■ (58) 



The averaged rate for a single fiavour-sensitive process is therefore [10 



3/(27r)2C(3)T3 ^-^^ • ^^^^ 

We have taken here for the numerical values g2 = 0.6 and gi = 0.4 at the time of 
Leptogenesis. Now with the order one approximation that all lepton doublets experience 
the same reaction rate, regardless of their momentum, we express 



such that it follows that 



36 

Bl = -^TIT = 1.7 X IQ-^T'^. (61) 



The difference when using a Fermi-Dirac instead of a Maxwell distribution in the above 
integral (164|) is, due to the factor |k| that favours contributions with large momentum, 
only about 6%. Using the same methods, we obtain for the weighted annihilation rate 
of Higgs particles 

Bl = 1.0 X lO-^T^ . (62) 
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As for the flavour sensitive rates, we proceed similarly and make our present estimates 



as follows: Ref. |46[ provides the number for the total production rate for sterile right- 
handed neutrinos. This gives an approximation for the flavour- sensitive scatterings that 
involve charged right-handed charged leptons. The main difference is that the charged 
right-handed leptons may also radiate U(l)y gauge bosons, an effect that is therefore 
neglected when adopting the number /i^/iS x 10~^T^ for the total reaction rate from 



Ref. j46|. The averaged rate for a single flavour-sensitive process then is 



'--wipm^"''''''~'^- '''' 

In the approximation of the same reaction rate for all leptons in phase space, we obtain 
Bfh^hSni = I ^|k|ri/.tM/, = / (^Ikle-I'^l/^ (64) 

and therefore 

Bf = ^T^T = 1.0 X lO-^Tl (65) 

In case a more accurate calculation of the lepton asymmetries is desired in the future, 
the estimates for the rates made in this Appendix, which should be of 0{1) accuracy, 
need to be improved. The main challenges are then to account for infrared effects and 
to abandon the approximations of averaging the reaction rates over phase space. 
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